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ABSTRACT 

We present measurements of the Type la supernova (SN) rate in galaxy 
clusters based on data from the Sloan Digital Sky Survey-II (SDSS-II) Super- 
nova Survey. The cluster SN la rate is determined from 9 SN events in a set 
of 71 C4 clusters at z < 0.17 and 27 SN events in 492 maxBCG clusters at 
0.1 < ^ < 0.3. We find values for the cluster SN la rate of (0.371°:}^!°:°}) SNur 

and (0.55i[3:}?+[3:[!?) SNur (SNux = lO-^^^^^^yr"^) in C4 and maxBCG 
clusters, respectively, where the quoted errors are statistical and systematic, re- 
spectively. The SN rate for early-type galaxies is found to be (0.31^o;}2iaoi) SNur 
h'^ and (0.49l°:ift°:°?) SNur h'^ in C4 and maxBCG clusters, respectively. The 
SN rate for the brightest cluster galaxies (BCG) is found to be {^MtlHlloZ) 
SNur h"^ and (0.36l!^J^l!^:|]J) SNur h'^ in C4 and maxBCG clusters, respectively 
The ratio of the SN la rate in cluster early-type galaxies to that of the SN la 
rate in field early-type galaxies is 1.94;j;^j+°;°^^ and Smtlfstomh for C4 and 
maxBCG clusters, respectively. The SN rate in galaxy clusters as a function of 
redshift, which probes the late time SN la delay distribution, shows only weak 
dependence on redshift. Combining our current measurements with previous 
measurements, we fit the cluster SN la rate data to a hnear function of redshift, 
and find tl = [(0.49l^:it)+ (0.9ll^|^) x z] SNuS /i^ A comparison of the radial 
distribution of SNe in cluster to field early-type galaxies shows possible evidence 
for an enhancement of the SN rate in the cores of cluster early-type galaxies. 
With an observation of at most 3 hostless, intra-cluster SNe la, we estimate the 
fraction of cluster SNe that are hostless to be (9.4; 



70. 
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Introduction 



The rate of Type la supernovae (SNe) in g alaxy clusters is an important field of study 
for a number of reasons. As discussed in, e. g., iGreggid (|2005[ ). supernova (SN) rate mea- 
surements are an observational probe of the progenitor systems, with the connection to 
progenitor models being made through inference of the distribution of delay times (DDTs) 
with respect to star formation. As galaxy clusters are generally composed of a high fraction 
of early-type galaxies that have old stellar populations, measurements of the SN la rate in 
clusters can in principle simplify the inference of the SN DDT. SNe in galaxy clusters are also 
a candidate source for metal enrichment of the intra-cluster medium (ICM). In particular, 
improved measurement of the rate of intra-cluster SNe would be significant for constrain- 
ing the relative importance of sources that may contribute to the cluster ICM enrichment 
(e. g. intra-cluster stars vs. galaxy outflow). 

The existing measurements of the cluster SN rate are few and are genera lly based on 



l ow-nu mbe r statisti c s. Est imates of the SN cluster rate were first presented by ICrane et al. 
( I1977I ) and iBarbonI (Il978l ). who c onsidered ^ 5 SNe discovered in the Coma cluster. The 



cluster S N rate was measured by iGal-Yam et al.l (120021 ) in clusters at ^ ~ 0.25 and z ~ 
0.9. The iGal-Yam et al.l (120021 ) results are based on a search for SNe in archival images of 
the Hubble Space Telescope (HST) and utilize one and two SNe, respectively. Subsequent 
to the HST SN search, a dedicated search for SNe in 161 Abel l clusters was undert aken 
by The Wise Observatory Optical Transient Search (WOOTS) dCal-Yam et~aDl2008h . A 
total of 6 SNe la di scovered by t h e WOO TS were used to determine the cluster SN la 
rate at z ~ 0.15 by ISharon et al.l (l2007b( ). A sample of 2-3 SNe la from the Supernova 
Legacy Survey (SNL S) have been used to determine the cluster SN la rate at z ~ 0.45 by 
Graham et al.l (12008 ). Fina lly, the sample of SNe in t he local jz < 0.0 4 ) univ erse presented 



by ICappellaro et al.l (Il999l ) have been reanalyzed by iMannucci et al.l ( l2008l ) to determine 
the cluster SN la rate with a samp le of 12.5 SN la (a fractional SN reflec ts uncertainty in 
typing; see ICappellaro et al.l ( 1l999l )). Additionally, IMannucci et al.l (120081 ) have placed the 
first constraints on the core-collapse (CC) SN rate in galaxy clusters based on a sample of 
7.5 CC SNe. A summary of the SN la rate results from the above cluster SN studies is given 
in Table [H A SN search in 15 massive, high-reds hift (0.5 < z < 0.9), X -ray selected clusters 
has been carried out on the HST, as described in ISharon et al.l ( l2007al ). and a measurement 



of th e SN la rate based on 6-14 SNe discovered by the program is forthcoming (ISharon et al. 
20091 ). A dedicated SN search, targeting ^ 60 X-ray selected clusters in the re dshift range 



0.1 < z < 0.2, is also being carried out on the Bok 2. 3m telescope on Kitt Peak (ISand et al. 
2008h 



As can be seen in Table H] the knowledge of the cluster SN la rate comprises 5 measure- 
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ments, based on a total of ~ 25 SNe. In this paper, we describe new measurements of the 
cluster SN la rate based on data from the Sloan Digital Sky Survey-II (SDSS-II) Supernova 
Survey (IFrieman et al.ll2008l ). The measurements are based on 35 SNe in the redshift range 
0.03 < z < 0.30, and therefore represent a significant statistical contribution to cluster SN 
la studies. In ^|2]we briefly describe the observations and SN search strategy of the SDSS-II 
Supernova Survey. In ^ we describe the galaxy cluster catalogs employed in this SN rate 
analysis. In ^we describe selection of the cluster SN sample from the SDSS-II Supernova 
Survey data. In §21 we describe necessary corrections to our SN la rate measurements. In 
^13 we present results on the cluster SN la rate, as well as limits on the cluster CC SN rate, 
and studies of the distribution of SNe with respect to their host galaxies. We summarize in 
^ Whenever necessary we assume a flat ACDM cosmology with Vt^ = 1 — = 0.3, and 
Hubble constant Hq = 70 km s~^ Mpc~^. 



2. SDSS-II Supernova Survey Observations 



Here we briefly describe aspects of the SDSS-II Supernova Survey most relevant to the 
present SN rate analysis. Much o f the materia l prese nted in this section is also relevant to 
the SN rate studies described in iDilday et al.l ( 120101 ) . and is discussed therein; we repeat 
the disc ussion here for the c onvenience of the reader. The survey is describ ed in more 



detail in iFrieman et al.l (|2008[ ) and the SN detection algorithms are described in ISako et al. 



(120081). Additional details of the survey observations an d the use of i n situ artificial SNe 
for determining SN detection effic iencies are discuss ed in iDilday et al.l (120081 ). A technical 



summary o f the SDSS i s give n bv lYork et a. 



be found in lHogg et al.l ( l200ll ): ISmith et al.l ( 



20021 



and quality assessment is describ ed by llvezic et al. 
described by iLupton et al.l (Il999[ ) 



(I2OOOI). Details of the su rvey calibration can 



Tucker et al.l ( l2006l ). the data processing 



( I2OO4J ). and the photometric pipehne is 



The SDSS-II Supernova Survey was carried out during the Fall (September- November) 
of 2005-2007, using the 2.5m telescope (iGunn et al.l l2006l ) a t Apache Point Obse rvatory 
(APO). Observations we re obtained in the SDSS ugriz filters ( jFukugita et al.lll996l ) with a 
wide- field CCD camera ( iGunn et al.lll998l ). operating in time-delay-and-integrate (TDI, or 
drift scan) m ode. The region o f the sk y covered by the SDSS-II Supernova Survey (designated 
stripe 82; see Istoughton et al.l ^00^ ) was bounded by -60° < ajaooo < 60°, and -1.258° < 
(^72000 < 1-258°. On average any given part of this ~ 300 square degree area was imaged 
once every 4 days during the survey operations. 



Difference images were produced in the SDSS gri filter bands by subtracting template 
images, constructed from previous observation of the survey region, from the survey images. 
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using an implementation of the methods described by lAlard &: LuptonI (jl998[ ). The difference 
images were sea rched for positive fluc tuations using the DoPHOT photometry and object de- 
tection package (ISchechter et al.lll993l ): typical limiting magnitudes (lOa above background) 
for the SDSS-II Supernova Survey were g ~ 21.8, r ~ 21.5, and i ~ 21.2. A combination of 
software cuts and human visual inspection was then used to identify promising SN candidates 
from the full set of transient detections. As a key component of prioritizing SN candidates 
for spectroscopic observation, the light curves for SN candidates were fit to models of Type 
la. Type Ib/c and Ty pe II SNe. This p rocedure is referred to as photometric-typing, and is 
described in detail by ISako et al.l ( 120081 ). 



Spectroscopic observations, for both SN type and redshift determination, were provided 
by a number of different telescopes. The spectra for the SNe utilized in the present SN 
rate analysis were provided by the Hobby-Eberly 9.2m at McDonald Observatory, the As- 
trophysical Research Consortium 3.5m at Apache Point Observatory, the William-Herschel 
4.2m, the Hiltner 2.4m at the MDM Observatory, the Subaru 8.2m on Mauna Kea, the 2.6m 
Nordic Optical Telescope and the 3.6m Italian Telescopio Nazionale Galileo at La Palma, 
the Mayall 4m telescope at Kitt Peak National Observatory, and the 3.5m ESO New Tech- 
nology Telescope (NTT) at L a Silla. Deta i ls of t he SDSS-II Supernova Survey spectroscopic 
data reductions are given by IZheng et al.l (|2008[ ). Comparison to high-quality SDSS galaxy 
spectra show that SN spectroscopic redshifts are accurate to ~ 0.0005 when galaxy emission 
features are used and ~ 0.005 when SN features are used. In either case the error on the 
spectroscopic SN redshifts are negligible for the SN rate studies considered here. 

While the difference imaging pipeline used during the SN search provides initial pho- 
tometric measurements, subsequent to the search more precise SN photometry is p r ovide d 
using a scene modeling photometry (SMP) technique developed by iHoltzman et al.l ( l2008l ). 
The final analysis of SN light curves and the selection cuts used to define the SN rate sample 
discussed in this paper are based on SMP. 



3. Galaxy Cluster Catalogs and Cluster Luminosity 



In studying the Type la SN rate in galaxy clusters we will work with two primary cluster 
catalogs; the C4 cluster catalog and the maxB CG cluster catalog. The C4 cluster finding 
algorithm and c atalog are disc u ssed in detail by iMiller et al.l ( l2005l ). The maxBCG catalog 
is presented by iKoester et al.l (l2007a{ ). and the cluster finding algorithm is described by 
Koester et al. ( 2007bl ). We briefly describe and summarize the content of these two cluster 
catalogs below. The redshift distributions for the clusters in these two catalogs that are in 
the SDSS-II Supernova Survey region are shown in Figured] 
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3.1. The C4 Cluster Catalog Description 



The C4 cluster catalog is based on the main spectroscopic sample of the SDSS and con- 
tains clusters in the redshift range 0.03 < z < 0.17. The C4 cluster catalog provides, among 
other quantities, cluster coordinates, the luminosity of each cluster in the SDSS r-band, and 
the number of galaxies identified as members of each cluster, A^^gais- The main spectroscopic 
sample of galaxies fr om SDSS is designed to be complete to a limiting magnitude of r ^ 17.8 
( IStrauss et al.ll2002[ ). For typical cluster galaxy luminosities, this implies that the identifica- 
tion of member galaxies, with L > 0.4 L^, is complete for clusters at z < 0.11. For clusters 
above this redshift limit a correction has to be made to the total cluster luminosity. The C4 
cluster identification algorithm works by searching for groups of objects that are tightly clus- 
tered in a 7-dimensional feature space, that includes spatial position, redshift, and observed 
colors. Note that there is no requirement that the colors for the galaxies be consistent with 
the colors of early-type galaxies, only that they be consistent with one another. The SDSS 
main galaxy sample cannot be 100% complete, as the SDSS fiber-spectrograph imposes a 
minimum angular separation of 55" (~ 100 kpc at z = 0.15) for objects targeted for 



the SDSS main spectroscopic galaxy sample (IStoughton et al.ll2002l ). A correction is applied 
to the C4 cluster luminosities to account for this by including galaxies as members of the 
cluster when they satisfy the magnitude requirement (r < 17.8) and have similar colors to 
the spectroscopically determin ed cluster members. T he published C4 catalog is based on the 
2nd data release of the SDSS ( lAbazajian et al.ll2004j ) and contains 748 clusters. The catalog 
used in this work is an extended version of th e C4 catalog that contains 1713 clusters and 
is based on the 5th data release of the SDSS ( lAdelman-McCarthy et al.l 120071 ) (B. Nichol, 
private communication). The subset of clusters that we will use in measuring the cluster SN 
rate are those that lie within stripe 82, and that have —40° < a;j2ooo < 50°. There are 71 
C4 clusters in this subset. 



3.2. The maxBCG Cluster Catalog Description 

The maxBCG catalog is based on SDSS photometric measurements, and the cluster 
identification algorithm relies on the tight relationship between color and redshift for lumi- 
nous red galaxies, which make up the majority of galaxy cluster composition. The maxBCG 
algorithm assigns a photometric redshift to each identified cluster that is derived by com- 
paring the cluster member galaxy colors to the expected colors for early-type galaxies, as a 
function of redshift. Comparison of the photometric redshifts for the maxBCG clusters to the 
spectroscopically measured redshift of the brightest cluster galaxy (BCG), when available, 
shows that the residuals for redshift (photometric - spectroscopic) are well described by a 
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Gaussian distribution with a mean of and a standard deviation of o" ~ 0.015. The observed 
{g — r) and (r — i) colors for maxBCG member galaxies, as a function of redshift, are shown 
in Figure |2J Linear functions were fit to to the (r — i) and {g — r) colors as a function of 
redshift, z, and these will be used below for applying k-corrections to the observations. For 
reference, the functions derived, valid for 0.1 < z < 0.3, are, 

(r-i) = 0.345 + 0.720 2 (1) 
(g-r) = 0.632 + 3.054 z. (2) 



The maxBCG catalog is restricted to the redshift interval 0.1 < z < 0.3, and is thus 
highly complementary to the C4 catalog. The lower-limit for the maxBCG cluster catalog 
is imposed because the fractional photometric redshift errors for redshifts less than z = 0.1 
are significant and have a large systematic effect on derived cluster properties. Furthermore, 
cluster catalogs based on spectroscopic redshifts for z < 0.1 are already available. The upper 
limit is imposed because, at z ~ 0.3, the "4000 A break" that is responsible for the uniformity 
in early-type galaxy colors moves into the region between the SDSS observer frame g and r 
band filters. Thus, the accuracy and precision for galaxy photometric redshifts is severely 
diminished ( iKoester et al.l l2007al ). A lower-limit for the luminosity of cluster members is 
imposed so that the definition of the composition of the clusters is consistent across the 
redshift range. The limit corresponds to an absolute magnitude of ^^■'^^^i ^ —20.25, (the 
(0.25)^ notation is explained in §3.3p and is such that the catalog is volume limited over 
the entire redshift range 0.1 < z < 0.3. Measures of the physical extent, r2005 and of the 
richness, A^'2oo and L20O) are provided for each cluster. r2oo is defined as the radius such 
that the mean density of early-typ e galaxies contained within is 200 times greater than 
the mean density of such galaxies (IKoester et al.l l2007al ): A^2oo(-^2oo) is defined as the the 
number (total luminosity) of early-type galaxies contained within r2oo- The public maxBCG 
cluster catalog contains 13,823 clusters with N20Q > 10. The subset of maxBCG clusters 
that we will use in measuring the cluster SN rate are those that lie within stripe 82, and that 
have —50° < aj2ooo < 60°. There are 492 maxBCG clusters in this subset. Additionally, 
this work makes use of a catalog of the maxBCG member galaxies (B. Koester, private 
communication) . 

In addition to the cluster richness and luminosity estimates that are pro vided with the 
maxBCG catalog, the raaxBCG clusters have been ext ensively studied by e.g. ISheldon et al. 
( l2009al ): IJohnston et al.l (120071 ): ISheldon et al.l ( l2009bl ). Of particular relevance to the study 
of the SN ra te in galaxy c l usters are the luminosity functions (LFs) of maxBCG clusters 
presented by iHansen et al.l ( l2009l ). and we will make extensive use of these below. 
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3.3. Luminosity Content of C4 and maxBCG Cluster Catalogs 

In ^we will present the SN rate in galaxy clusters per unit luminosity. In this section 
we present a comprehensive discussion of the luminosity content for the C4 and maxBCG 
cluster catalogs. In what follows we use the notation ^^^m to denote the observer frame 
magnitude, m, for an object that has been k-corrected to a redshift z. All luminosities 
presented here are for galaxies with L > 0.4 L^, {L^ is a characteristic luminosity for cluster 
members; see Equation [3]). Correcting for the faint end of the luminosity distribution is 
discussed below. 

The C4 catalog provides total cluster luminosities in SDSS r-band. The total uncorrected 
r-band luminosity for the galaxies identified as cluster members, for clusters considered in 
this study, is 2.02 x 10^'^ Lq /i"^, where Lq is the luminosity of the sun, and h is the value of 
the Hubble constant in units of 100 km s~^ Mpc~^. As discussed in §3.H these luminosities 
have to be corrected to account for incompleteness of the spectroscopic sample due to fiber- 
collisions. The total r-band luminosity after making this correction is also provided by the 
C4 catalog and has the value 4.08 x 10^^ Lq h~'^. A mentioned above, the luminosities 
for clusters at z > 0.11 also have to be corrected to account for the fact that some cluster 
galaxies with L > 0.4L^ will have observed magnitudes fainter than the completeness limit 
of the main SDSS galaxy sample (r ^ 17.8). Of the 71 C4 clusters we are considering, 12 
are at z > 0.11. Applying a co rrection to the lumin osities of these 12 clusters, assuming the 



cluster luminosity functions of iHansen et al.l ( 120091 ). results in a total r-band luminosity for 



C4 clusters in this study of 4.12 x 10^^ L© h~'^. 

The maxBCG catalog includes the summed luminositi es for membe r galax ies in 



(0.25)„ 



and (0-25) However, the maxBCG cluster LFs presented by [Hansen et al.l ( 20091 ) represent a 
more complete study of the luminosity content of maxBCG clusters, and we will use these as 
the definitive measure of the total maxBCG cluster luminosities. In particular, a background 
subtraction has been performed that reduces inaccuracies in the cluster luminosities due to 
interloping fo reground and backgr ound galaxies that may be counted as maxBCG cluster 



members. In [Hansen et al.l ( l2009l ). luminosity functions (LFs) are presented for maxBCG 
satellites, as a function of the richness measure, iV2oo- The LFs account for both red and 
blue cluster galaxies, but do not include the contribution to the luminosity from the BCG. 
The LFs are assumed to take the form of a Schechter function, 

mdL = 4>* (AV e-^/^-^ (3) 



where is a characteristic luminosity for cluster members, and (J)* is a normalization con- 
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stant with units of inverse volume. The 3 parameters of the LFs, (j)*, L^,a, are each expressed 
as functions of A'^200) with the general functional form of A x A^^^qq. The values of A (normal- 
ization) and /9 (exponent) for these 3 LF parameters are given in Table [2J To use the LFs 
to compute the total luminosity we exploit the identity 

E L = {L{N2oo)) N200 (4) 

where EL denotes the summed luminosity for the cluster, and the average luminosity, (L), 
is given by, 

r„°?r dL L (j){L) 

Equation m gives the total C^-^^)^ luminosity in a cluster, as a function of A''2oo- The total lumi- 
nosity in clusters in the SN survey region is then the sum over iV2oo of {L(N2oo)) NsNiN2oo), 
where Nsn{N2oo) is the number of clusters in the survey region for which the number of 
member galaxies is A'"2oo- Using this formalism we find that the total luminosity of maxBCG 
clusters in the survey region is SLf^^ = 1.096 x 10^^ Lq h~^. 



3.3.1. Correcting Cluster Luminosities for Faint Galaxies 

The luminosities quoted above for the C4 and maxBCG cluster catalogs include only 
galaxies with L > 0.4 L*, which is a conventional way of characterizing cluster luminosities. 
In measuring the cluster SN rate we do not wish to exclude SNe occurring in faint galaxies, 
and so it is necessary to estimate the contribution to the total cluster luminosities from 
galaxies with L < 0.4 L*. The total cluster luminosity can be estimated as L = k where 
denotes the luminosity for galaxies with L > 0.4 and the correction factor, k, is given 

by 

IoIl, ^ ^(^) 

For a typical maxBCG cluster with A^2oo = 20 the power-law exponent of the luminosity 
function is a = —0.59 (Table |2]). For a cluster with A^2oo = 35 the value for the exponent is 
a = —0.68. The corresponding correction factors are k = 1.21 and k, = 1.25. We will assume 
that the faint end of the luminosity function is a characteristic property of galaxy clusters 
and that the power-law behavior of the LFs for maxBCG clusters is appropriate for the C4 
clusters also. 
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Table 1. Cluster Rate Measurements 



Reference 


Redshift 


Mean 


Lookback Time NsNe 


SN la Rate 




Range 


Redshift 


[Gyr] 


[SNuB h^] 



This work (C4) 


0.03 


- 0.17 


0.084 


1.11 


9 





46 


This work (maxBCG) 


0.10 


- 0.30 


0.225 


2.69 


25 





68 


Mannucci et al. f2008) 


- 


0.04 


0.020 


0.28 


12.5 





57 


Sharon et al. f2007b) 


0.06 


- 0.19 


0.150 


1.89 


6 





73 


Gal- Yam et al. f2002) 




0.25 


0.250 


2.94 


1 





80 


Graham et al. f2008) 




0.45 


0.450 


4.67 


3 





63 


Gal- Yam et al. f2002) 




0.9 


0.900 


7.30 


2 


1 


63 



+0.21+0.01 
-0.15-0.01 

+0.17+0.02 
-0.14-0.02 



-0.16 

+0.45 
-0.29 

+ 1.84 
-0.65 

+ 1.04 
-0.33 

+2.16 
-1.06 



'^For this work, the quoted errors are statistical and systematic, respectively. For pre- 
vious measurements, the total error is quoted; see the corresponding references. 



Table 2. maxBCG Luminosity Function Parameters 



LF Parameter 


Normalization 


Exponent 




(A) 






8.0 Mpc"^ 


-0.20 




0.8 X 10^° Lq 


0.15 


a 


-0.28 


0.25 



Note. — Par ameters refer t o Sch echter func- 
tions derived by [Hansen et al.l (120091 ) . 



3.3.2. k-corrections 



In order to minimize uncertainties due to k-corrections, the maxBCG magnitudes, and 
the corresponding luminosities, have been k-corrected to the median redshift of the clusters, 
z = 0.25. For comparison of the SN rate results based on the maxBCG clusters with 
previous cluster SN rate measurements, it will be necessary to k-correct the luminosities 
into more standard filters. In order to determine the appropriate k-corrections we selected a 
set of galaxies from the SDSS galaxy catalog that satisfy the color vs. redshift relations for 
maxBCG members discussed above. The corresponding set of galaxies contains ~ 675, 000 
members. We then matched these galaxies with their corresponding records in the photoz 
database provided by the SDSS catalog arc hive server (CAS), which provides k-corrections 
based on the work of Blanton et al. J2003ah . We thereby determine an average k-correction 



as a function of redshift, appropriate to the early- type galaxies that make up the bulk of the 
maxBCG catalog. The k-corrections so derived, in SDSS r and i bands, are shown in Figure 
ini For reference, the functions employed for k-corrections as a function of redshift, z, valid 
for 0.1 < z < 0.3, are 

Krr{z) = 9.17 X 10"^ + 1.48 z (7) 
KJz) = 3.95x10-3 + 1.012. (8) 



The functions above were derived by binning the k-corrections in redshift and taking the 
error on the mean (root mean square divided by square root of the number of entries) as 
the uncertainty. With this formalism the error on the fitted parameters is < 1% and is thus 
negligible in comparison to the statistical error on our SN la rate measurements. 

To transform the maxBCG cluster luminosities from (O-^^)^ to r, we note the following 
identity: 

nir = m^^-'^^ - Krr{z = 0.25) (9) 
= mf^^ + {r -i){z = 0.25) - Krriz = 0.25) (10) 



We k-correct the maxBCG cluster luminosities from (O-^^)? to r because at a redshift z = 0.25 



rest-fr ame r maps closely to observer-frame i, and because the analysis of [Hansen et al. 



(120091 ) that we take as the definitive measurement of the maxBCG cluster luminosities is 
performed for ("-^s)^ and not ('^■'^^^r. Using the expressions given above for (r — i) and Krr, 
evaluated at z = 0.25, we derive the transformation = '•^■^^^mj + 0.063. For comparison 
to other cluster SN rate measurements, it is necessary to express the cluster luminosities in 
units of Lq. The conversion to solar luminosities is given by 
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L 



L 



^ -^g-O.4((m.-mo.25,)-(M5-M0 ')) 



'11^ 



where Mq is the absolute magnitude of the sun. To compute the absolute magnitude of the 
sun in an arbitrary filter, we use a solar spectrum obtained from the CALSPEC0 database, 
hosted by the Space Telescope Science Institute, and compute synthetic magnitudes us- 
ing a custom piece of software written for this task. We thereby derive values for the 
absolute magnitud e of th e sun of Mq = 4.64 (in agreement with the value published by 
Blanton fc Roweid (120071 )). and Mq ' = 4.67. With these values. Equation [U] becomes 
{L/LQ)r = 0.918 (L/L0)o.25j, giving the r-band luminosity in maxBCG clusters, for galaxies 
with L > 0.4 L„ as SL^ = 1.006 x 10^^ Lq. 



3.4. Uncertainty on the Cluster Luminosities 

Any uncertainty on the luminosity of the clusters described above will propagate to 
an uncertainty on the derived cluster SN la rates. Sources of uncertainty on the cluster 
luminosities might include, e.g., uncertainties in the background subtraction used in defining 
the maxBCG cluster luminosities, or uncertainty in the correction for fiber collisions in 
defining the C4 clusters. The cluster catalogs described above do not include uncertainties 
on the cluster luminosities, and it is beyond the scope of this analysis to determine these 
uncertainties. However, as an estimate the order of magnitude of the cluster luminosity 



uncertainties, we note that iMenanteau et al.l ( 120101 ) give uncertainties on the luminosities of 



individual clusters, for a cluster catalog that was constructed in a similar way to the maxBCG 
cluster catalog. These uncertainties are ~ 10%, and include both statistical and systematic 
effects. As discussed above, the SN la rate analysis includes 71 and 492 clusters from the 
C4 and maxBCG catalogs, respectively. If the uncertainties on the individual clusters are 
~ 10% and are independent, then the uncertainty on the total cluster luminosities will be 
< 1%. We conclude that uncertainty on the cluster luminosities is likely to be negligible 
in comparison to the statistical uncertainty on the cluster SN la rate measurement, but we 
acknowledge that any additional uncertainties on the cluster luminosities, or correlations 
between the uncertainties on the luminosity of the individual clusters, would impact the 
precision of our cluster SN la rate measurements. 



littp://www. stsci.edu/hst/observatory/cdbs/calspec. html 
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Fig. 1. — Redshift distributions for the C4 {z < 0.17) and maxBCG (0.1 < z < 0.3) clusters 
in the SDSS-II SN Survey region. There are 71 clusters from the C4 catalog and 492 clusters 
from the maxBCG catalog. 

4. The Cluster SN Sample 

4.1. Type la Supernovae 

To define the SN sample that is associated with galaxy clusters it is necessary to define 
selection criteria on both the SN light-curve properties, and the SN spatial correlation with 
clusters. As was done in ou r study of the low redshift SN rate from the first season of the 



SDSS-II Supernova Survey (jPilday et al.ll2008l ). the SNe in the rate sample are required to 



meet certain sel ection criteria on their light-curve fits to the MLCS2k2 light curve model 



(jJha et al.l 120071 ). For the present analysis the selection criteria are, 



At least 1 observation at < —2 days relative to peak in the supernova rest-frame. 

At least 1 observation at > tiate days relative to peak in the supernova rest-frame. 
For maxBCG cluster SN candidates tiatc = 10 days, and for C4 cluster SN candidates 
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Fig. 2. — Observed colors of maxBCG galaxies vs. redshift. The points represent the mean 
at the corresponding redshift, and the error bars represent the error on the mean. The solid 
lines show the best fitting linear function. 

tiate = 5 days. The peak of the SN light curve is defined as the epoch of maximum 
luminosity in the rest-frame i?-band. 

• Fit probability > 10~^. 

The fit probability is computed assuming that the value of the minimum from the 
MLCS2k2 light curve fit follows a distribution with Ai'p — 4 degrees of freedom, 
with iVp the number of photometric observations. All SNc light curves arc fit for 
epoch of maximum light, luminosity parameter (A) and extinction parameter (Ay). 
When the SN redshift is known from spectroscopic measurement, the SN light curve is 
additionally fit for distance modulus. In the case of photometric SNe, the additional 
fit parameter is the redshift of the SN. 



The requirement of a late time photometric measurement is used primarily to reduce misiden- 
tification of photometric SN la candidates. For the low redshifts of the C4 galaxy clusters 
the majority of SNe have spectroscopic confirmation of type, and so a requirement of tiate > 5 
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Fig. 3. — K-corrections vs. redshift for early- type galaxies. The points represent the mean 
at the corresponding redshift, and the error bars represent the error on the mean. The solid 
lines show the best fitting linear function. 

days increases the size of the SN sample, without introducing uncertainty from misidentifi- 
cation of photometric SNe. For the maxBCG clusters, the number of photometric cluster SN 
candidates is significant and a stricter requirement, tiate > 10 must be used. Additionally, 
the SNe that will be used for determining the SN rate are required to satisfy a set of selection 
criteria, based on fits of their search photometry to models of Type la. Type II and Type 
Ib/c light curve i nodels. These addit ional selection criteria are identical to those discussed 
in more detail in iDilday et al.l ( 120101 ). Briefly, the SN candidates must be discovered in at 
least 3 epochs of the SN search, must fit the Type la model well in relation to the Type II 
and Type Ib/c models, an d must not have an overly broad light c urve shape, as typified by 
the peculiar SN la 2005gj (lAldering et al.ll2006l : iPrieto et al.ll2007l ). The combination of cuts 
on light-curve sampling and shape effectively reject AGN from the cluster SN la sample. 

A SN candidate is defined to be associated with a galaxy cluster if it satisfies 



SN is within 1 Mpc h ^ projected distance of the center of a cluster. 
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• SN redshift is consistent with the cluster redshift. 

Projected distance here refers to the orthogonal distance from the SN to the center of the 
galaxy cluster, assuming the redshift of the cluster. The choice of 1 Mpc h^^ projected 
distance to define association with a cluster is chosen largely for consistency with previous 
cluster SN la rate measurements. The definition of redshift consistency depends on whether 
the SN and/or cluster redshifts are determined photometrically or spectroscopically. For 
the C4 clusters, the cluster redshifts are always precisely determined with a spectroscopic 
measurement, whereas the maxBCG cluster redshifts are determined photometrically. Con- 
sistency between the SN redshift, 2;^, and cluster redshift, is defined in the following 
way: 

• Spectroscopic SN redshift and spectroscopic cluster redshift 
\zs - Zc\ < 0.015 

• Spectroscopic SN redshift and photometric cluster redshift 
1^, - Zc\ < 0.025 

• Photometric SN redshift and spectroscopic cluster redshift 
1^, - Zc\ < 2.5 V(0.01)2 + §^2 

• Photometric SN redshift and photometric cluster redshift 
1^, - Zc\ < 2.5 V(0-015)2 + 5z^ 

where Szg is the error on the SN photometric redshift. In the case that both the SN and the 
cluster have spectroscopically determined redshifts, the allowable spread in galaxy redshifts 
is taken as 0.015. This number was determined from the catalog of maxBCG cluster mem- 
ber galaxies. For each cluster that has at least 4 member galaxies with spectroscopically 
determined redshifts, we computed the root mean square (RMS) of the difference between 
the member spectroscopic redshifts and the redshift assigned to the cluster. The mean of the 
resulting distribution of RMS values is ^ 0.015. A similar value is obtained if the difference 
between the cluster redshift and the member galaxy spectroscopic redshifts are instead fit 
with Gaussian distributions, and the mean of the distribution of RMS values for each Gaus- 
sian is computed. We note that for the 9 C4 cluster SNe that have spectroscopic redshifts for 
both the SN and the cluster, the deviation in redshift is < 0.002 and this cut could be made 
more strict without affecting the SN rate result. In the case that the SN has a spectroscopi- 
cally measured redshift and the cluster has a photometrically determined redshift (maxBCG 
clusters only) the tolerance on redshift consistency (0.025) is determined by adding in quadra- 
ture the RMS spread of maxBCG member galaxy redshifts mentioned above (0.015) and the 
typical accuracy of the maxBCG cluster photometric redshifts (?» 0.015 — 0.020). 
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In the case that the SN redshift is determined photometrically, the cuts we apply cor- 
respond to a tolerance of 2.5 a^, where represents the error on the SN redshift, added in 
quadrature to a tolerance on the cluster redshift. This is a rather loose tolerance on redshift 
consistency, and we discuss this issue further in §3 There are 9 SNe la in C4 clusters and 27 
SNe la from maxBCG clusters that satisfy the selection criteria, and these are listed in Table 
[3] (C4 clusters) and Table H] (maxBCG clusters), respectively. We note that SNe 12979 and 
18375 occur in a cluster that is a member of both the C4 and maxBCG catalogs, and thus 
there are 34 distinct SNe la that pass the selection criteria. We note that SNe 16280, 18047 
and 18362 are each associated with two distinct clusters, and that SNe 14279 and 16215 are 
associated with the same cluster. In the cases where a SN is associated with two clusters, 
we will count it only once in the SN rate calculation, assigning it to the nearest of the two 
clusters. 



Corrections to the SN la Rate Measurements 



5.1. Search Efficiency 



The method for determinin g the SN discovery efficiency is based on the same Monte 



Carlo (MC) studies discussed in iDilday et al.l (|2010[ ). However there is an important mod- 
ification to be considered for SNe in galaxy clusters. It is a well established result that 



Sullivan et al. 


(2006 


); 


Jha et al. 


(2007) 



SNe (e. g. 

sumed distribution for the luminosit y parameter o f the M LCS2k2 light curve model. A, for 
the entire SN sample, as employed in iDilday et al.l (120101 ) . is not appropriate for considering 
SNe in galaxy clusters. To determine the SN discovery efficiency for SNe in galaxy clusters, 
we generated a set of MC SNe with a A distribution similar to the observed distribution 
for cluster SNe. The A distribution assumed for cluster SNe la is shown in Figure ID The 
fiducial A distribution has mean (A) = 0.033, with an RMS of 0.286, and the cluster SNe 
A distribution has mean (A) = 0.081, with an RMS of 0.246. Assuming a larger mean 
value for A does not significantly affect the SN selection efficiency at the low redshifts of 
the C4 clusters, and has a 2% effect at the redshifts of the maxBCG clusters. The SN 
selection efficiency at the redshifts of the C4 clusters is approximately constant for each 
observing season, with the values 0.77, 0.73 and 0.72 for the 2005, 2006 and 2007 observing 
seasons, respect ively. The re d shift dependent efficiencies used for the maxBCG cluster SNe 
are discussed in IDilday et al.l ( 12010[ ). Briefiy, the efficiencies are wel l described by a func tion, 
en/(lie^o)/flK with eo ^ 0.7, Zo ^ 0.35 and s, ^ 0.05. As in IPildav et all toO^ and 



Dilday et al.l ( l2010l ) the efficiency at low-redshift is significantly less than 1 due mainly to 
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Fig. 4. — Assumed distribution for the MLCS2k2 A parameter for the entire SN sample 
(black/sohd) and for the intrinsically fainter cluster SNe (gray/dashed). 



the requirements of early and late time observing epochs. 



5.2. Contamination Prom Redshift Uncertainties 



In §U we have defined a set of criteria on the consistency of the SN and cluster redshifts, 
which, in the case of photometric redshifts, correspond roughly to a requirement that the 
redshifts match to within ^ 2.5cr. We chose to apply a loose cut on the photometric redshift 
consistency because cluster SNe are rare by their nature and we wish to maximize the size of 
the sample, and because we assert that the requirement of a spatial correspondence with a 
cluster is already a strong filter on contamination from chance projections. Here we quantify 
further the expected contamination rate from chance projections, in the presence of finite 
error on the SN redshift measurement. We note that for the C4 cluster SN sample, all the 
SNe have spectroscopically measured redshifts (see Table [3]), such that the contamination 
rate due to using SN photometric redshifts is manifestly 0, and we focus exclusively on 
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contamination in the maxBCG SN sample. 

To estimate the expected contamination rate we employ the following procedure. For 
each cluster in the maxBCG catalog we compute the expected contamination, 

dzv(z) dz'^^eM \ , ^ ), (12) 

where //(z) is the SN number density, per unit redshift, z^x is the cluster redshift, is the 
assumed error on the SN photometric redshift, and n^- defines our cut on redshif t consistency. 



he e rror on the photometric redshift, as a function of redshift, is shown in iDilday et al. 
(|20101) to be ~ 0.2 x z^'^ . To compute ?7(-2), we use the power-law rate model from 



Dildav et all mm [ry = 2.6 x 10^^ (1 + z)^-^ SNe Mpc"^ h^o), along with the SDSS-II 
Supernova Survey SN efficiency function. The solid angle used in computing ri{z) is the solid 
angle corresponding to a 1 Mpc separation from a cluster at redshift z^. However, we 
assume that no SNe will be spatially associated with a cluster and have a redshift within 
±0.015 of the cluster unless it is, in fact, associated with the cluster. Therefore ri{z) is 
assumed to be for —0.015 < z — Zd < +0.015. The choice of 0.015 as the spread of 
maxBCG cluster redshifts is discussed in §4.11 For each cluster in the maxBCG catalog, we 
also compute the number of expected cluster SNe la, N'^, assuming a cluster SN la rate pe r 



unit r-band lum i nosity that is broadly consistent with those given by ISharon et al.l (l2007bl ): 



Mannucci et al.l (|2008[ ). The observed number of SNe that would be included in the sample. 



A^*^, is then related to the true number, A^"^, by 

]\fO ]\T 

^ = erfK/^2)) + ^^ (13) 

where erf(.) is the error function, and Ai'contam is the total expected contamination, Sj nl^^^^^^. 
To clarify our approach, we note that the right-hand side of Eqn. (fT3|) can be computed, 
with reasonable assumptions about the value of the cluster SN la rate. We can then com- 
pare this to the observed value of the left-hand side of Eqn. (fT3l) . and thereby infer the 
contamination rate of the SN la sample due to finite redshift uncertainties. For cuts on 
redshift consistency of = 1.5 and n^- = 2.5, A^contam/-^"^ has the values 0.151 (rd/SNur)"-^ 
and 0.276 (rd/SNur)^^, where is the assumed value of the cluster SN rate, and SNur = 
(10^° Uq)-\100 yt)~\ 

We first note that using values for no- ~ 1.5 — 2.5, we consistently arrive at A^^ ~ 5. By 
way of example, if we choose = 2.5 and Vd = 0.6 SNur, then we find /N'^ = 1.44 (r.h.s. 
of Eqn. ( !T3|) ). There are 7 SNe la with photometrically determined redshifts that are within 
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±2.5(7 of the cluster redshift (Table [3]), and so this implies a value of A^"^ = A^*^/1.44 ^ 5 
(l.h.s. of Eqn. f lT3|) ). Therefore, we apply a bias correction to our maxBCG SN sample by 
assuming that the number of photometric SNe in the sample is 5, rather than the 7 listed in 
Table HI To estimate the systematic error on the size of the maxBCG SN sample, we var y 
the assumed power-law SN rate model according to the errors given in iDilday et al.l ( l2008l ). 
and vary the assumed cluster rate as = (0.6 ± 0.2) SNur, which results in an estimated 
systematic error of ~ (—2 — 1-3)% on the maxBCG cluster SN rate measurements, due to 
employing photometric redshifts for a subset of the cluster SNe. 



5.3. Cluster Incompleteness and Contamination 

Here we discuss the effect of incompleteness and contamination in the galaxy cluster 
catalogs on our cluster SN rate measurements, and on our comparisons of the cluster SN rate 
to the rate in early type galaxies in the field. In what follows we assume that the SN rate in 
both galaxy clusters and the field is proportional to the host luminosity, with proportionality 
constants etc and a/ for clusters and the field, respectively. We n ote that the C4 c luster 



catalog has ~ 10% incompleteness and an ^ 5% contamination rate (IMiller et al.ll2005l ). and 



that the maxBCG ca talog has ~ 10% incompleteness and an ^ 10% contamination rate 



flKoester et al.ll2007b[ ) 



Cluster Incompleteness 

With the ansatz that the cluster SN la rate is proportional to the stellar luminosity, 
independent of cluster richness, cluster incompleteness does not affect the derived value of 
the SN rate. That is, although real clusters may fail to be included in the cluster catalog, 
the cluster SN sample will be reduced in direct proportion to the luminosity of the missing 
clusters. Cluster incompleteness, however, will affect the value of the ratio of the rates 
for cluster vs. field early-type galaxies, as SNe occurring in clusters may be erroneously 
associated with the field SN sample. 

In the presence of cluster incompleteness, the inferred ratio of the SN la rate in clusters, 
Tc, to the SN rate in field early-type galaxies, Vf, will be. 



Nc-Ns-\ f Lf + Ls- 



Tf \L,-Ls-J \Nf + Ns- 



(14) 



^(1 + (1 - aJaf){Ls-/L,)iLjLf)), (15) 
af 

where N (L) represents the true number of SNe (luminosity) in a given subset, and the c. 
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/, and 6— subscripts denote clusters, the field, and the portion of the clusters missed due to 
incompleteness, respectively. The quantity Ls-/Lc is the cluster incompleteness. To estimate 
the quantity L^/Lj, we compute a cluster luminosity density by dividing the summed r-band 
luminosity from the maxBCG cluster catalog (k-corrected from as described in 



for clusters in the SN survey region, by the corresponding volume. For redshift s z > 0.2 this 
derived cluster luminosity density has the value ~ 0.104 x 1 0^ L^d Mpc~^ h . In 



Di 



day et al. 



( 120081 ) we showed that, using the luminosity functions of iBlanton et al.l ( l2003bl ). and an 
objective classification scheme for early-type galaxies, the "'^r-band luminosity density, for 
early type galaxies is 0.994 x 10^ Lq Mpc~^ h. Applying a k-correction as described in 
§3.3[ this corresponds to a r-band luminosity density, for early- type galaxies, of 1.24 x 
10^ Lq Mpc~^ h. Therefore, we conclude that the ratio of the total r-band light in cluster 
galaxies vs. field early- type galaxies, LrJL f, is ^ 0.091. The fraction of cluster luminosity 



contained in early-type galaxies is ~ 80% ( iHansen et al.ll2009[ ). and thus the ratio of the r- 
band light in cluster early-type galaxies vs. field early-type galaxies, (Lc/L/)red, is ~ 0.073. 
In turn, this results in an estimated — (1 — 2)% correction factor, for values of the ratio of 
the SN rates, Oc/a/, of ~ 2 — 3. 

Cluster Contamination 

Contamination of the cluster catalogs will affect both the value of the SN rate in clusters, 
and the ratio of the rate in clusters to the rate in the field. 



The inferred cluster SN rate will be. 



Lc + Ls+ 

ac (1 + {af/otc 



1) Ls+/L, 



(16) 
(17) 



where N and L are defined as in equation [TH and 5+ denotes field quantities that were erro- 
neously identified with clusters. The quantity L^^/Lc represents the cluster contamination 
rate. The estimated correction factor on the SN la cluster rate measurement due to cluster 
contamination is then +3% for C4 clusters, and ~ -|-(5 — 7)% for maxBCG clusters. 

The inferred ratio of the SN rate will be, 

r, ^ f N, + Ns+ \ / Lf-Ls+ \ 

^ —{l + {af/ac-l)Ls+/L^), 
af 

with estimated correction factors on the cluster vs. field SN rates of ~ +3% for C4 clusters, 
and ~ +(5 — 7)% for maxBCG clusters. 
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5.4. Summary of SN la Rate Corrections 

To summarize our consideration of necessary corrections to our cluster SN la rate mea- 
surements, we have considered the effects of the lesser average intrinsic luminosity for SNe 
la in early-type galaxies as compared to a representative galaxy sample, contamination of 
the cluster SN la sample due to the use of SN photometric redshifts, and incompleteness 
(clusters are not identified by the cluster finding algorithms) and contamination (objects that 
are not clusters are identified as clusters) of the galaxy cluster catalogs. The lesser average 
intrinsic luminosity of SNe la in early-type galaxies does not affect the selection efficiency 
for the relatively low-redshift C4 SN sample, and has an ^ 2% effect on the maxBCG SN 
sample. This factor is accounted for through the efficiency function, e{z), that is mentioned 
in ^ We estimate that ^ 2 of the 7 SNe with photometric redshifts in the maxBCG SN 
la sample are not in fact associated with the cluster, and thus we will apply a correction 
of /^^ ^ —2/27 = —7.4% to the maxBCG cluster rate. Cluster contamination results in 
correction factors of ~ +3% for the C4 cluster SN la rate, and ~ +(5 — 7)% for 
the maxBCG cluster SN la rate. Cluster contamination also results in correction factors 
on the ratio of the SN la rate in early-type cluster galaxies to early-type field galaxies of 
^ +3% for the C4 clusters, and /'^ ~ +(5 — 7)% for the maxBCG clusters. Additionally, 
cluster incompleteness results in an ^ —(1 — 2)% correction factor on the ratio of the 
SN la rate in early-type cluster gal axies to early- t ype fi eld galaxies for both the C4 and the 



maxBCG clusters. We note that in iDilday et al.l (120101 ) it is shown that the contamination 
of the SDSS-II SN Survey photometric SN la sample by non-la SNe is not more than ^ 3%. 
Photometric SNe la make up ~ 50% of the maxBCG cluster SN la sample, and so the overall 
contamination of the cluster SN la sample is < 2%. Contamination of the SN la sample due 
to lack of a spectroscopic identification is likely to be correlated with the contamination due 
to the use of photometric redshifts discussed above, and so the size of the effect is expected 
to be less than this. 



6. Systematic Errors 

Here we discuss possible sources of systematic error on our cluster SN la rate measure- 
ments. In the previous section we derived correction factors to account for incompleteness 
and contamination of the galaxy cluster catalogs. These correction factors will each have 
an uncertainty due to sample variance in the galaxy clusters in the SN survey region. We 
evaluate the uncertainty by assuming that each galaxy cluster in the sample will suffer from 
incompleteness or contamination with a probability given be the mean incompleteness and 
contamination rates given above. The uncertainty on the correction factor is then derived 
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from the variance of the corresponding binomial distribution, with the number of events equal 
to the number of galaxy clusters considered. For the C4 clusters (71 clusters), this results in 
a ±3.2% uncertainty on the correction for cluster incompleteness and a ±2.5% uncertainty 
on the correction for cluster contamination. For the maxBCG clusters (492 clusters), this 
results in a —1.2% - ±1.0% uncertainty o n the correction for both cluster incompleteness 
and cluster contamination. iDilday et al.l (120101 ) discuss the systematic uncertainty on the 
SN la rate for the full SDSS-II SN Survey sample, and show that assuming a value for the 
mean extinction in V^-band, (Ay) = 0.45, as opposed to the fiducial value of (Ay) = 0.35, 
has a large systematic effect on the derived SN la rates for SNe la at z > 0.2. For the present 
sample of SNe la in galaxy clusters, assuming a larger than fiducial value for the mean dust 
extinction is not appropriate, and so we do not include this systematic effect on the cluster 
SN la rate measurements. As noted in §3.41 any uncertainty on the cluster luminosities will 
also contribute to the total uncertainty on the cluster SN la rate. 



7. SDSS SN Results 

7.1. The C4 cluster rate 

As discussed above, there are 9 type-la SNe in C4 clusters that satisfy the SN selection 
criteria. The total r-band luminosity in the SN survey region, after correcting for the faint 
end of the LF, is 4.99 x 10^^ Lq /i"^, and the mean redshift of the C4 clusters is {z) = 0.0786. 
The observing time for the SDSS-Il Supernova Survey was 89, 90, and 90 days, for the 2005, 
2006, and 2007 observing seasons, respectively. The SN selection efficiency is approximately 
constant over the range of the C4 clusters, with the values 0.77, 0.73, and 0.72 for the 2005, 
2006, and 2007 observing seasons, respectively. Using these values, the SN rate in C4 clusters 
is 

r = 1^^^ = 0.37^°:}^lHi SNur (20) 

where i denotes each observing season, T is the survey observation time, e is the SN selection 
efficiency, N is the number of SNe la, SL^ is the total cluster luminosity in the r band, and 
SNur = (10^° Lq)~^(100 yr)~^. The factor = 1 + = 1.03, where is the correction due 
to contamination of the cluster catalog, discussed in ^ The errors quoted are the 1-sigma 
statistical and systematic errors, respectively. We have used the mean survey efficiency, 
determined from Monte Carlo simulations, as both the positions of the cluste rs and the 



Monte Carlo SNe are effectively uniform r a ndoin samplings of the survey area. iJorgensen 



( 119971 ) (clusters) and IPadmanabhan et al.l ( 120041 ) (field) give the average stellar mass to 



luminosity ratio (r-band) for early-type galaxies as ~ 3, which is in good agreement with 
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the M/L conversion employed by ISharon et al.l (l2007bl ). Using this assumption, the Type 
la SN rate per unit luminosity quoted above is equal to 0.123;!;q;q4q1q;qq3 SNuM h'^, with 

Using the same average conversion factor from r to -B as 



SNuM = (10 ^ ° M^)- HlOO yr)- 



Sharon et all f l2007bh . this corresponds to 0.46t°:2it0;°} SNu5 h'^. 



7.2. The maxBCG cluster rate 



As discussed in §11 there are 27 SNe in maxBCG clusters from SDSS-II Supernova 
Survey that satisfy the selection criteria. As the SDSS-II Supernova Survey SN discovery 
efficiency is not well approximated as constant over the redshift range of the maxBCG 
catalog, in determining the SN rate in maxBCG clusters we use the more formal definition 
of the SN rate per unit luminosity, r^. 



TL = (21) 
eTSL 



with 



-^min 

where is the earth frame observation time, ULi^z) is the total cluster luminosity as a 
function of redshift, and ei^z) is the SN discovery efficiency. The factors g^^^ and g'^ represent 
the corrections due to use of SN photometric redshifts and contamination of the cluster 
catalog, respectively i^^. These factors have the values g^^ = 1 + = 0.926 and g'^ = 
1 + = 1.06. For the maxBCG catalog z^am is fixed to 0.1. In Figure |5] we show the value 
of the maxBCG cluster SN rate as a function of ^max- This figure shows that the derived 
SN rate is not strongly sensitive to the upper limit on the cluster sample. If the upper limit 
is chosen as the upper limit of the maxBCG catalog, z = 0.3, then we have N = 27 and 
eTHL = 4.85 x 10^^ yr Lq h~'^, including the correction for the faint end of the LF. The 
derived value of the SN rate is thus tl = O.^btofitofi SNur (= 0.68l|]:JIl|]:|J^ SNu5 
= 0.18;o:o36lo:oo4 SNuM h^). Figure [6] shows the cluster SN rate as a function of the limit 
on the projected distance from the center of the cluster, in units of Mpc h^^. The fraction of 
red, early-type galaxies in clusters is larger at small separations from the cluster center, and 
in this sense the SN rate at smaller separations is a more reliable probe of the component 
of the SN rate that originates from an old stellar population. Furthermore, the extent of a 
cluster is not an unambiguously defined quantity, and SNe at smaller separations are more 
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robustly associated with the cluster. Figure [6] shows that the derived SN rate is not strongly 
dependent on the limiting projected radius that we use to define cluster membership. Figure 
[7] shows the SN rate as a function of the lower limit on the cluster richness measure, A^2oo- 

The SDSS-II Supernova Survey cluster rate results, along with the previous measure- 
ments listed in Table [H in units of SNu-B /i^, are shown in Figure [H A fit of the data to a 
linear model of the cluster SN la rate as a function of redshift, r = A + Bz, gives best fit 
values ofA = 0.49 ±0.14 SNuE and B = 0.8llJ]|^ {x'^/NDF = 0.88/5). A fit of the data 
to a constant model of the cluster SN la rate as a function of redshift, r = A, gives a best 
fit value ofA = O.QltoH SNu5 {x^/NDF = 1.9/6). 



7.3. Cluster SN Rate vs. Field SN Rate 



In iDilday et al.l (120081 ). the SN rate in low-redshift early-type galaxies was estimated 
from the first year data of t he SDSS-II Superno va Survey as ~ 0.17]'^q;q4 SNur h"^. One SN 
in the sample described in iDilday et al.l (120081 ) is a cluster SN, and hence the early-type 



field rate is ^ 0.16l|J:[!^ SNur h'^ Using the s ame c olor cut [u — r = 2.4) to differentiate 



early and late type galaxies as IDilday et al.l (120081 ). we find that 6/9 of the C4 cluster 



SNe reside in early-type galaxies, while 20/27 maxBCG cluster SNe reside in early-type 
galaxies (18/25 afte r bias correction to account for SN photometric redshifts). According to 
Hansen et al.l (120091 ). the fraction of light in early- type galaxies for the maxBCG clusters is 
~ 80%. Therefore the cluster SN rate for early-type galaxies is = O-^^^-o.il-om SNur h"^ 
and = 0.491o;}^1o;q^ SNur h'^ for C4 and maxBCG clusters, respectively. As discussed 
in ^ the ratio of the cluster early-type to field early-type rates must be multiplied by a 
factor = 1 — = 0.985 to account for incompleteness of the galaxy cluster catalogs. 
The C4 and maxBCG cluster SN la rates are seen to be larger by factors of 1.94lQ;9}lQ;Qf5 
and 3.02l};g3lo;Q4g compared to the field SN rate for early-type galaxies, respectively. This 
is broad agreement with th e enhancement of the Type la SN rate in early-type galaxies in 
galaxy clusters reported by lMannucci et al.l ( 120081 ). 



7.4. The SN la Rate in Cluster BCGs 

Of the 32 SNe in our cluster sample (after accounting for duplicates and for contami- 
nation due to use of SN photometric redshifts), 2 are in the brightest cluster galaxy (BCG) 
of the cluster. These are SNe 12979 and 13905. Another, SN 18890 appears to be in the 
BCG, although it has a somewhat ambiguous host. All 3 of these SNe are in C4 clusters. 



-27- 




Fig. 5. — maxBCG cluster SN rate as a function of upper limit on the redshift range. The 
dashed lines represent the 1-sigma upper and lower limit of the SN rate. 

while 1 is in a maxBCG cluster. The cluster BCGs contain ^ 6% of the total r-band cluster 
luminosity, and so the SN la rate, in the cluster BCGs, is ri = 2.04;;*;};^^i'^q;q4 SNur h'^ for the 
C4 clusters and = 0.36^°:®^+°;°} SNur h"^ for the maxBCG clusters. 



7.5. Hostless SNe 



For the SNe la in our cluster sample, all but 2 are visually associated with a host galaxy. 
The 2 SNe la with no host evident in the images are SN 13073 and SN 19001. Additionally, 
one SN, SN 1782, lies at a distance ~ 25 kpc from its presumed host galaxy, and might 
reasonably be considered to be hostless. Assuming SN 1782 is in fact a hostless SN, we 
can constrain the fraction of hostless SNe in gal axy clusters to be ~ 9.4%, with 1-sigma 



confidence interval (4.3 — 17.7)%. For comparison, iGal-Yam et al.l (120031 ) reported 2 hostless 
SNe from a sample of 7 for an estimated fraction of hostle ss SNe of ~ 30%. We note. 



however, that the classification of these 2 SNe as hostless by iGal-Yam et al.l (120031 ) relied 
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Fig. 6. — maxBCG cluster SN rate as a function of limit on the projected distance. The 
dashed lines represent the 1-sigma upper and lower limit of the SN rate. 



in part, on the observation of velocity offsets between the SNe and its potential host of 
^ 750 — 2000 km/s. In principle a more thorough analysis of SN and galaxy spectra for 
SNe from the SDSS-II Supernova S urvey may reveal additional hostless SNe. The clusters 
considered by lGal-Yam et al.l ( 120031 ) were generally more massive than the clusters considered 
here, and an additionally possibility is that the rate of hostless SNe la in galaxy clusters is 
larger for more massive clusters. 



7.6. Radial Distribution of Cluster SNe 



Forster fc Schawinskil (120081 ) have studied the radial distribution of SNe in early-type 



galaxies. The sample of SNe includes an unspecified number of SNe from the SDSS-II 
Supernova Survey. Here we consider the radial distribution for SNe for field and cluster 
ellipticals. 



To construct the sample of early-type galaxies, we employ the following procedure. We 
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Fig. 7. — maxBCG cluster SN rate as a function of lower limit on the cluster richness 
measure, A^2oo- The dashed lines represent the 1-sigma upper and lower limit of the SN rate. 



first assign a host galaxy to each SNe by locating the nearest host from the SDSS database, 
in units of the isophotal radiu s of the host galaxy in r-band. The algorithm is described 
in detail in iDilday et al.l (|2008[ ). Early-type galaxies are defined as those host galaxies that 
satisfy the following criteria: 



M - r > 2.4 
r < 21.5 
Ar < 0.05 



where u and r are the SDSS model magnitudes and Ar is the error on the r-b and magnitude. 
It is a well established result that the {u—r) color for SDSS galaxies is bimodal (IStrateva et al. 



20011) . with early-type galaxies generally having {u — r) > 2.2. As discussed in 



Dilday et al. 



(I2OO8I ). st udies of the observe d u — r distribution for galaxies from the photometric redshift 
catalog of lOyaizu et al.l (120081 ) suggest that u—r = 2.4 may provide a more robust separation. 
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The requirement that r < 21.5 is imposed as the separation of stars and galaxies is fairly 
robust to this limit. The requirement that Ar < 0.05 is imposed to remove outlying, poorly 
measured galaxies. The distributions of the distance of each SN from its host galaxy, for 
field and cluster early-type hosts, is shown in Figure M For field early-type galaxies, a fit 
of the data to a Sersic model of the luminosity distribution, dN/dr = A p e~'^^ ^, where p 
denotes the distance of the SN in units of the deVaucouleurs radius of the host galaxy, gives 
a value for A of 0.20±0.08, which is consistent with a deVaucouleurs profile (A = 0.25). Since 
the distribution of light in early-type galax ies is known to follow a deV aucouleurs profile, 
the result of the fit confirms the results of iForster fc Schawinskil (120081 ) that the SN rate 
in field ellipticals is well represented by a constant rate per unit luminosity. The radial 
distribution for SNe in cluster early-type galaxies shows an an enhancement at small radial 
separations in comparison to a deVaucouleurs profile. An enhancement of the SN la rate 
in regions that have und ergone rec e nt sta r formation has been reported by several authors 
(see, e.g., the review by iMannuccil (120091 )). and our observed enhancement at small radial 
separations is possible evidence for a component of the cluster SN la rate that tracks residual 
star formation activity in cluster early-type galaxies. Such an enhancement is in qualitative 
agreement with the larg er SN la rate in cluster ellipticals, compared with field ellipticals, 
mentioned above and by lMannucci et al.l (l2008l ). as well as with recen t independent e viden ce 



for some rec ent star-formation activ ity in early-type galaxies (e.g, iKaviraj et al.l ( 120071 )). 
Additionally, iDella Valle et al.l (|2005[ ) and iGraham et al.l (|2010[ ) report an enhancement of 
the SN la rate in radio-loud early-type galaxies, and attribute this to mergers that not only 
power the radio emission, but also provide a young stellar population that can account for the 
enhanced SN la rate. While we do not consider the radio properties of the e arly-type galaxies 



i n our sam ple, it is possib l e that our results are related to the results of iDella Valle et al. 



(120051 ) and iGraham et al.l (120101 ). While the best fit radial profile for cluster early type 
galaxies does not match a deVaucouleurs profile, a KS test on the radial distributions of 
SNe la in cluster and field ellipticals results in an 30% pr obability tha t the data are 
drawn from the same underlying distribution. As discussed in iDilday et al.l ( 120081 ). all SN 
selection efficiency calculations are based ultimately on artificial SNe la inserted directly 
into the survey search imaging data. Analysis of these artificial SNe la from the three 
observing seasons of the SDSS-II Supernova Survey does not show evidence for significant 
loss of efficiency near the cores of galaxies. While the artificial SNe were inserted into random 
galaxies that on the average may be less luminous in the cores than typical cluster galaxies, 
any additional inefficiency would only increase the observed enhancement of the SN la rate 
in the cores of cluster early-type galaxies. 
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Fig. 8. — Cluster SN la rate vs. redshift. The dashed hne shows the best fit to a hnear 
model of the SN rate as a function of redshift. 

8. Conclusions 

We have presented measurements of the Type la SN rate in galaxy clusters over the 
redshift range 0.03 < z < 0.30. These measurements are based on 34 SNe la (32 after 
applying a bias correction) from the SDSS-II Supernova Survey and represent a significant 
statistical contribution to the study of the Type la SN rate in galaxy clusters. Our results 
on the Type la SN rate are consistent with previously published measurements, both in the 
local universe, and at redshift ^ 0.15 — 0.25. In turn, the local and low-redshift SN cluster 
rates are consistent with the SN la rate at z ^ 0.45 and at z ^ 0.9. The current cluster 
SN la rate measurements do not show evidence for a rapid increase in the SN rate, as a 
function of lookback time, as has been well-established for the volumetric SN la rate. It 
should be emphasized, however, that the existing cluster SN rate measurements are based 
on small samples and the measurements do not rule out a redshift dependence to the cluster 
SN rate. A fit of the cluster SN la rate measurements to a linear dependence on redshift 
results in a best fit slope of (O.QI^q g^) SNuS /i^ per unit redshift. We find a ratio of the SN 
la rate in cluster early-type galaxies to that of the SN la rate in field early-type galaxies of 
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Fig. 9. — Radial distribution of SNe in early- type host galaxies. 



1.941^:;^}+°:°^^ and 3.02t}i^^°:°^^^ for C4 and maxBCG clusters, respectively. 

We find that at most 3 of the SNe la in our sample arc hostless, intra-cluster SNe, 
which is significantly less than the 30% hostless fraction estimated from previous cluster SN 
studies. 

We have presented the first study of the radial distribution of SNe la in cluster early-type 
galaxies. The radial distribution for SNe in cluster early-type galaxies shows an enhancement 
at small radial separations in comparison to the radial distribution in field early-type galaxies, 
which are well described by a deVaucouleurs profile. This enhancement could be attributable 
to residual star formation in cluster early-type galaxies, which could explain the higher SN 
la rate observed in cluster early-type galaxies in comparison to early-type galaxies in the 
field. 
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Note. — SDSS Id denotes internal candidate designation. The horizontal Unes delineate SNe discovered during the 3 distinct observing seasons of the SDSS-II Superijcwa 
Survey. 00 

''"la" refe rs to spectroscopica lly identified SN; see IZheng et al. I ||2008| '). ' 'la-photo+z" refers to photometrically identified SN with a spectroscopically measured host galaxy 
redshift; see lPildav et al.l 1I2OIOI) . 

°SN does not satisfy the color- typing criteria. See iDildav et al. I ||2010|') . 

''Value has been k-corrected to the rest frame of the galaxy; see i|3.3.2l 

'^Observer frame magnitudes for the SN host galaxy. 



Table 4. SNe in maxBCG Clusters. 
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Note. — SDSS Id denotes internal candidate designation. The horizontal lines delineate SNe discovered during the 3 distinct observing seasons of the SDSS-II Supernova Survey. 

^ "la" and "la?" refer to spectroscopically ide ntified SNe; seelZheng et al. "la-photo+z" and "la-photo" refer to photometrically identified SNe with and without a spectroscopically 

measured host galaxy redshift, respectively; see iDildav et al. 1 ||2010|'I . 



alue has been k-corrected to the rest frame of the galaxy; see i|3.3.2l 

rrors are shown for photometric SNe la. For spectroscoipcally measured redshifts, the error is negligible and is not listed ( ^2i 
observer frame magnitudes for the SN host galaxy. 
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